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ERF  IS  A  POTENTIAL  ERK-MODULATED  TUMOR  SUPPRESSOR  IN  PROSTATE 
CANCER 


INTRODUCTION 


Metastatic  prostate  cancer  is  the  second  most  frequent  cause  of  cancer  death  among 
American  men  and  remains  incurable.  Half  of  such  cancers  contain  a  gene  fusion,  placing  a 
slightly  truncated  ERG  transcription  factor  downstream  of  the  androgen-regulated  promoter  of 
TMPRSS2.  This  leads  to  significant  upregulation  of  ERG,  which  drives  tumorigenesis.  ERG 
belongs  to  the  ETS  family  of  transcription  factors,  which  contain  the  highly  conserved  ETS 
DNA-binding  domain.  There  are  differences  in  the  biology  of  fusion-positive  tumors  versus 
those  lacking  the  ERG  fusion  (fusion-negative).  However,  there  is  no  targeted  clinical  strategy 
for  treating  either  fusion-positive  or  fusion-negative  cancers.  This  is  partly  due  to  our  limited 
understanding  of  how  exactly  ERG  promotes  oncogenesis,  as  well  as  a  lack  of  awareness  of 
complimentary  pathways  in  fusion-negative  tumors.  A  recent  multi-institutional  effort 
sequencing  metastatic  prostate  cancers  yielded  ERF  as  1  of  only  19  genes  achieving  the  highest 
statistical  criteria  for  containing  cancer-relevant  mutations,  with  several  frameshift  mutations  as 
well  as  point  mutations  of  key  conserved  residues  noted  in  the  critical  ETS  domain,  presumably 
corresponding  to  loss-of-function  mutations.  ERF  has  known  tumor  suppressor  ability  and 
belongs  to  the  same  transcription  factor  family  as  ERG,  contains  a  virtually  identical  DNA- 
binding  ETS  domain,  but  is  a  strong  transcriptional  repressor  unlike  the  transcriptionally 
activating  ERG.  Most  interestingly,  all  tumors  with  ERF  mutations  were  fusion-negative.  We 
hypothesize  that  in  fusion-negative  tumors,  ERF  functions  as  a  tumor  suppressor  by  occupying 
ETS  binding  sites  and  inhibiting  gene  expression,  consistent  with  its  repressor  function.  On  the 
other  hand,  in  fusion-positive  tumors,  loss  of  ERF  function  confers  no  discernible  growth 
advantage,  since  ERF  is  already  outcompeted  at  ETS  binding  sites  by  the  upregulated  ERG 
oncogene,  which  in  turn  activates  gene  transcription. 


KEYWORDS 


prostate  cancer,  tumor  suppressor,  TMPRSS2-ERG 
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ACCOMPLISHMENTS 


Goals  and  Accomplishments 

I.  Specific  Aim  1 :  To  investigate  how  ERF  and  ERG  may  compete  and  have  opposing 
effects  on  cancer  cell  proliferation 

A.  Major  Task  1 :  Evaluate  role  of  ERF  in  negatively  regulating  ERG-dependent  genes 
Nothing  to  report  (please  see  2016  annual  report  for  previous  data). 


B.  Major  Task  2:  Evaluate  role  of  ERF  in  competing  for  DNA  binding  with  ERG 

Last  year,  I  described  the  development  and  analysis  of  a  VCaP-based  system  to  study 
competition  between  ERG  and  ERF  binding  for  chromatin.  These  are  two  ETS  transcription 
factors  that  have  a  similar  DNA-binding  motif  but  have  intrinsically  positive  or  negative 
transcriptional  activity,  respectively.  The  VCaP  system  was  ideal  for  initial  studies,  since  it  has 
high  levels  of  endogenous  ERG.  However,  this  left  open  the  possibility  that  the  competition  that 
we  observed  was  a  cell-line  specific  phenomenon,  and  not  generalizable  to  other  prostate  cancers 
or  normal  prostate  tissue. 

To  this  end,  and  given  the  role  that  ERF  also  plays  in  benign  prostate  cells,  I  performed  ChlP- 
seq  to  confirm  that  ERF  binds  ETS  sites  in  normal  prostate  organoids.  It  is  important  to  note  that 
ERG  is  not  endogenously  expressed  in  normal  prostate,  unlike  half  of  prostate  cancers  (Fig.  la). 
De  novo  motif  analysis  of  the  ERF  ChIP-seq  peaks  identified  the  canonical  GGAA  ETS  DNA 
motif  as  the  primary  ERF  binding  site  (Fig.  le).  Remarkably,  transient  overexpression  of  ERG 
led  to  a  significant  decrease  of  ERF  chromatin  occupancy  (Fig.  la-e),  consistent  with  the 
competition  for  binding  also  seen  previously  in  VCaP  cancer  cells.  Furthermore,  the  ChIP-seq 
signal  intensity  of  the  smaller  number  of  ERF  peaks  observed  in  the  FLAG-ERG  state  was 
increased  in  the  FLAG-alone  state  in  almost  all  cases  (Fig.  lb).  Sh3bp5  and  Rph3al  are  example 
genes  demonstrating  the  ERG-mediated  inhibition  of  ERF  chromatin  occupancy  (Fig.  Id). 
Given  the  similarity  of  the  ETS  domains  between  ERF  and  ERG,  as  well  as  the  ability  of  ETS 
factors  to  displace  one  another  from  canonical  GGAA  motifs,  we  postulated  that  the  opposing 
effects  of  ERF  and  ERG  could  be  explained  by  competition  for  ETS  binding  sites.  I  explored  the 
overlap  between  androgen  receptor  and  ERF  sites  in  normal  prostate  organoids  by  androgen 
receptor  ChIP-seq  and  found  28%  of  the  ERF  sites  overlap  with  androgen  receptor  binding  sites 
(Fig.  le).  This  was  lower  than  observed  in  the  previously  analyzed  ERG-positive  VCaP  tumor 
cells,  consistent  with  differences  in  androgen  receptor  chromatin  occupancy  between  normal 
prostate  and  cancer.  Collectively,  these  data  suggest  a  model  whereby  the  relative  levels  of  ERF 
and  ERG  compete  for  ETS  binding  sites. 
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Figure  1  -  a)  Mouse  normal  prostate  organoids  were  infected  with  a  TetOn  doxycycline-inducible  Flag-ERG  or  empty  vector 
(Flag-alone)  and  treated  with  or  without  doxycycline.  b)  ERF  ChIP-seq  in  normal  prostate  organoids  infected  with  either  Flag- 
alone  or  Flag-ERG  lentivirus,  both  treated  with  doxycycline  (n  =  2  biological  replicates:  Rl,  R2),  and  analysed  with  heat  maps,  c) 
Comparison  of  ERF  ChIP-seq  peak  numbers,  n  =  2  biological  replicates  (Rl,  R2).  d)  ERF  ChIP-seq  signals  for  example  genes 
Rph3al  and  Sh3bp5  loci,  normalized  by  read  depths,  e)  ERF  ChIP-seq  in  normal  prostate  organoids  infected  with  either  a  Flag- 
ERG  lentivirus  or  Flag  alone,  both  treated  with  doxycycline  (n  =  2  biological  replicates:  Rl,  R2)  analysed  by  (1)  peak  overlap,  (2) 
motif  analysis,  and  (3)  overlap  with  ERG  and  androgen  receptor  ChIP-seq. 
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C.  Major  Task  3:  Evaluate  for  in  vitro  role  of  shERF 


The  near  complete  lack  of  ERF  mutations  in  tumors  with  TMPRSS2-ERG  fusions  led  us  to 
investigate  whether  ERF  loss  recapitulates  the  phenotype  of  ERG  gain,  keeping  in  mind  that 
ERG  is  not  expressed  in  benign  prostate  epithelium.  Last  year,  I  described  the  development  and 
analysis  of  normal  murine  prostate  organoids  with  ERF  knockdown.  However,  ERG  only 
displays  potent  oncogenicity  in  a  mouse  Pten- A  background.  To  this  end,  and  to  determine 
whether  ERF  loss  yields  a  phenotype  similar  to  ERG  gain,  I  now  infected  prostate  organoids 
derived  from  Pten ~h  mice  with  short  hairpin  RNA  (shRNA)  targeting  mouse  Erf  (shErf  m).  The 
Pten-A  shErfm  organoids  acquired  morphological  characteristics  of  ERG  overexpression:  they 
formed  luminal  structures  lacking  basal  cells  (Fig.  2a),  and  profiling  by  RNA  sequencing  (RNA- 
seq)  and  qPCR  demonstrated  profound  loss  of  expression  of  most  basal  signature  genes  (Fig. 
2b-e).  Consistent  with  these  changes,  RNA-seq  profiling  revealed  significant  enrichment  for 
genes  whose  expression  is  upregulated  by  androgen  in  human  prostate  cancer  cells  (Fig.  2c  - 
right  panel),  which  will  be  explored  further  in  Specific  Aim  3.  Collectively,  this  supports  the 
hypothesis  that  ERF  loss  leads  to  a  pro-luminal  shift,  consistent  with  the  changes  observed  as 
normal  prostate  tissue  transitions  to  adenocarcinoma. 
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Figure  2  -  a)  Immunohistochemistry  of  Pten  '-  mouse  prostate  organoids  infected  with  non-targeting  shRNA  (shNT)  or  targeting 
ERF  (shErf_m).  b)  Basal  signature  applied  to  Pten+/+  organoid  RNA-seq  (n  =  2  biological  replicates),  c)  RNA-seq  analysis  of 
organoids  derived  from  Pten+/+  and  Pten1-  mouse  prostates  infected  with  non-targeting  shNT  or  shErf  m;  n  =  2  biological 
replicates,  d)  RT-qPCR  analysis  of  the  Pten+/+  organoids;  n  =  2  biological  replicates,  e)  Pten+/+  organoid  RNA-seq  (n  =  2 
biological  replicates)  interrogated  by  GSEA  for  expression  signature  of  Witte  basal  prostate  cancer 
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Specific  Aim  2:  To  evaluate  for  in  vivo  tumor  suppressor  function  of  ERF 
A.  Major  Task  1 :  Evaluate  role  of  ERF  as  in  vivo  tumor  suppressor 

The  Pten~h  organoids  infected  with  shErf_m  were  indeed  able  to  form  tumors  when  grafted  back 
into  mice,  recapitulating  the  phenotype  of  ERG  overexpression  and  confirming  that  ERF  is  a 
bona  fide  tumor  suppressor  (Fig.  3a).  Next,  to  assess  the  impact  of  ERF  in  human  prostate 
cancers,  we  interrogated  the  mRNA  expression  profiles  of  the  primary  and  metastatic  human 
prostate  cohorts  TCGA-333  and  SU2C-150  respectively.  In  agreement  with  the  prior  functional 
studies,  ERF  mRNA  levels  are  inversely  correlated  with  two  luminal  transcriptional  activity 
signatures  both  in  normal  human  prostate  and  in  all  primary  tumor  subtypes  (Fig.  3b).  This 
reciprocal  association  is  also  observed  in  metastatic  cancers  if  the  analysis  is  limited  to  tumors 
without  amplification  of  or  mutations  in  the  AR  gene,  which  encodes  the  androgen  receptor  (Fig. 
3c).  The  fact  that  reduced  ERF  expression  enhances  luminal  transcriptional  output  even  in  the 
absence  of  mutation  or  deletion,  raises  the  possibility  that  ERF  may  have  a  broader  role  in 
prostate  oncogenesis. 
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Figure  3  -  a)  Graft  tumour  volumes  of  GEMM-derived  prostate  organoids.  Data  are  median  ±  interquartile  range;  n  =  10  tumours 
per  condition,  b)  Expression  profiles  of  the  TCGA-333  primary  prostate  cancer  cohort  (n  =  333  patients)  were  interrogated  for 
correlation  between  the  ERF  mRNA  level  and  two  androgen  transcriptional  activity  signatures.  P  values  were  calculated  by  the 
Spearman  correlation  test,  c)  The  same  analysis  as  b  was  applied  to  the  SU2C-150  (n  =  150  patients)  metastatic  castration-resistant 
prostate  cancer  cohort3  (mCRPC). 
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B.  Major  Task  2:  Evaluation  of  role  of  point  mutants  in  ERF  tumor 


We  queried  all  published  prostate  cancer  genome  cohorts  (n  =  930  patients)  and  complied  a  list  of 
all  ERF  mutations  (Fig.  4a),  present  in  1-3%  of  patients.  The  mutations  include  the  specific 
K401fs  and  G299fs  loss-of-function  truncations  (initially  defined  in  craniosynostosis  families), 
as  well  as  similar  missense  mutations  in  the  ETS  domain.  Mapping  of  these  ETS  missense 
mutations  onto  the  known  crystal  structure  of  ERG  revealed  that  the  altered  residues  are  located 
within  conserved  helices  (Fig.  4b),  with  four  of  the  five  conserved  mutations  predicted  to  be 
destabilizing  by  in  silico  analysis.  Expression  of  complementary  DNAs  (cDNAs)  containing  the 
ETS  missense  or  truncating  mutations  in  LNCaP  prostate  cancer  cells  led  to  reduced  ERF  steady- 
state  levels  relative  to  wild-type  ERF  despite  robust  messenger  RNA  (mRNA)  expression  (Fig. 
4c,  d),  in  agreement  with  prior  craniosynostosis  studies.  Collectively,  this  pointed  to  ERF  point 
mutations  in  prostate  cancer  as  being  destabilizing. 
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Figure  4  -  a)  ERF  point  mutations  identified  in  all  prostate  cancer  cohorts.  *Mutations  shared  with  patients  having 
craniosynostosis.  Green:  ETS  missense  mutations  of  residues  conserved  with  ERG.  Red:  frameshift  (fs)  /  splice-site  (sp)  mutations 
upstream  of  repressor  domain,  b)  ERG  ETS  domain  (Protein  Data  Bank  accession  number  4IRI)  illustrating  sequence  conservation 
with  ERF.  ERF  missense  mutations  indicated  by  sticks,  c)  Expression  of  ERF  mutant  proteins  in  LNCaP  cells,  d)  ERF  mutant 
protein  densitometry  from  immunoblot  in  c  compared  with  mRNA  RT-qPCR. 
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C.  Major  Task  3:  Evaluation  of  ability  of  MEK/ERK  inhibitor  to  augment  tumor  suppressor 
activity  of  ERF 

Nothing  to  report 


Specific  Aim  3:  To  interrogate  and  map  how  ERF  may  regulate  androgen-mediate 
transcription  via  its  novel  interaction  with  AR 

A.  Major  Task  1 :  Determination  of  role  of  ERF  in  AR  signaling 

The  inverse  correlation  between  ERF  mRNA  level  and  the  androgen  receptor  signature  is  also 
observed  in  the  half  of  prostate  cancers  that  possess  the  TMPRSS2-ERG  fusion  (noted 
previously  in  Fig.  3b,  c),  suggesting  that  ERF  may  have  an  androgen  receptor  repressive  function 
in  this  subtype  as  well.  To  investigate  this  possibility,  the  expression  of  ERF  (Fig.  5a)  or  ERG 
was  separately  inhibited  via  shRNA  in  the  ERG-fusion-positive  VCaP  cell  line  and  the  androgen 
transcriptome  analyzed,  looking  for  genes  with  a  two-fold  change  in  expression  by  RNA-seq 
upon  InM  dihydrotestosterone  (DHT)  addition,  FDR<0.05.  Consistent  with  prior  literature, 
inhibition  of  ERG  expression  (‘ERG-low’)  resulted  in  a  contracted  androgen  transcriptome 
compared  with  the  wild-type  state  (‘ERG-high’)  (Fig.  5b  -  Left  panel).  Conversely,  ERF 
inhibition  (‘shERF_2’)  increased  the  change  in  expression  of  androgen  receptor  target  genes  and 
doubled  the  size  of  the  androgen  transcriptome  (Fig.  5b  -  Left  panel).  DHT  treatment  had  no 
effect  on  ERF  mRNA  levels  (Fig.  5b  -  Mid  panel)  and  doxycycline  alone  had  no  effect  on  gene 
expression  (Fig.  5b  -  Right  panel).  ERF  knockdown  had  no  effect  on  ERG  protein  level,  on 
androgen  receptor  protein  levels,  or  on  androgen  receptor  subcellular  localization  (Fig.  5a). 
Consistent  with  data  observed  in  Fig.  2c  and  Fig.  3b-c,  ERF  knockdown  in  VCaP  cells  led  to 
enrichment  of  androgen-signaling  genes,  both  noted  by  GSEA  interrogation  of  the  above  RNA- 
seq  experiment  (Fig.  5c)  or  by  validation  RT-qPCR  experiments  (Fig.  5d).  Collectively,  ERF 
serves  as  an  endogenous  repressor  of  AR  signaling,  both  in  normal  prostate  and  in  prostate 
cancer  cells. 
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Figure  5  -  a)  Nuclear/cytoplasmic  fractionation  of  VCaP  cells  infected  with  a  non-targeting  shRNA  (shNT)  or  shRNA  targeting 
ERF  (shERF_2).  b)  Left  -  Androgen-regulated  genes  (at  least  a  two-fold  change,  FDR  <  0.05  by  RNA-seq  with  1  nM  DHT  for 
16h)  in  VCaP  cells  infected  with  a  doxycycline-inducible  shRNA  targeting  ERG  (ERG-low)  or  targeting  ERF  (shERF_2)  and 
analysed  by  heat  map;  n  =  3  biological  replicates.  Center  -  Analysis  of  same  RNA-seq  data  evaluating  the  effect  of 
dihydrotestosterone  (DHT)  on  ERF  expression.  Right  -  Analysis  of  RNA-seq  data  to  evaluate  the  effect  of  doxycycline  addition 
alone,  c)  Interrogation  of  RNA-seq  in  shERF_2  VCaP  cells  (n  =  3  biological  replicates)  by  GSEA  for  Nelson  androgen  up 
expression  signature,  d)  RT-qPCR  of  shERF_2-infected  VCaP  cells  treated  ±  DHT.  Data  are  mean  ±  s.e.m.;  n  =  3  biological 
replicates. 
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B.  Major  Task  2:  Confirmation  of  AR-ERF  interaction 

C.  Major  Task  3:  Significance  of  AR-ERF  interaction 


Nothing  to  report 


Training  Opportunities  and  Professional  Development 

I  have  been  undergoing  an  intensive  training  program  for  a  prospective  prostate  cancer 
physician-scientist.  I  continue  to  attend  bioinformatic  seminars  necessary  for  the  RNA-seq  and 
ChIP-seq  analysis  in  Aims  1  and  2.  I  attend  oncology  grand  rounds  weekly.  I  also  attend  the 
interdisciplinary  bimonthly  prostate  cancer  meeting  organized  by  our  Genitourinary  Service,  and 
attend  the  monthly  inter-institutional  prostate  cancer  meeting  held  amongst  the  prostate  cancer- 
focused  labs  in  New  York  City.  Regarding  conferences,  I  attended  the  Prostate  Cancer 
Foundation  meeting.  Dr.  Howard  Scher  continues  to  provide  clinical  mentorship  regarding  my 
work  with  patients,  and  Dr.  Charles  Sawyers  continues  to  be  my  mentor  with  respect  to  this 
research  project. 


Dissemination  of  Results 

The  work  to  date  was  published  (please  see  ‘Products’  section). 


Next  Reporting  Period 

Aiming  to  make  progress  on  the  clinical  information  (e.g.  prognosis  and  prediction  to  therapy 
response)  associated  with  patients  whose  tumors  possess  ERF  mutations. 
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IMPACT 


On  principal  discipline 

The  findings  over  the  reporting  period  have  furthered  our  understanding  of  tumor  genetics 
responsible  for  prostate  cancer  progression.  Specifically,  this  project  has  investigated  the  poorly 
understood  nature  of  the  ERG  fusion-negative  tumors  and  further  clarified  the  nature  of  the  ERG 
fusion-positive  tumors  as  well.  The  identification  of  an  ETS  factor,  ERF,  that  serves  as  a  tumor 
suppressor  in  fusion-negative  tumors,  but  is  outcompeted  in  ERG  fusion-positive  prostate  cancer, 
is  a  significant  contribution  to  our  understanding  of  the  tumor  genetics  of  prostate  cancer. 

The  results  from  this  proposed  research  project  could  also  have  a  major  impact  on  prostate 
cancer  care.  Existing  treatments  for  advanced  prostate  cancer  increase  survival  only  on  the  order 
of  months,  and  can  be  associated  with  significant  morbidities.  Half  of  all  such  cancers  lack  ERG 
gene  fusions,  and  at  present,  we  lack  a  unifying  understanding  of  the  key  regulators  in  such 
tumors.  If  the  project  hypothesis  is  correct,  then  we  may  be  able  to  identify  a  novel  treatment 
strategy  for  patients  possessing  fusion-negative  tumors.  This  approach  could  lead  to  increased 
survival  and  decreased  morbidity  for  men  suffering  from  this  disease.  In  addition  to  the  clinical 
applications  above,  this  project  could  shed  light  on  how  tumors  with  and  without  ERG  mutations 
may  actually  be  similar.  This,  in  turn,  would  contribute  to  our  understanding  of  the  50%  who  do 
have  the  ERG  mutation,  and  assist  in  developing  a  targeted  strategy  for  them  as  well. 


On  other  disciplines 
Nothing  to  report 


On  technology  transfer 
Nothing  to  report 


Beyond  science  and  technology 
Nothing  to  report 


CHANGES/PROBEEMS 


Nothing  to  report 
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PRODUCTS 


Publications,  conference  papers,  and  presentations 
Publications: 

Bose  R,  Karthaus,  WR,  Armenia  J,  Abida  W,  Iaquinta  PJ,  Zhang  Z,  Wongvipat  J,  Wasmuth  EV, 
Shah  N,  Sullivan  PS,  Doran  MG,  Wang  P,  Patruno  A,  Zhao  Y,  International  SU2C/PCF  Prostate 
Cancer  Dream  Team,  Zheng  D,  Schultz  N,  Sawyers  CL.  ERF  mutations  reveal  a  balance  of  ETS 
factors  controlling  prostate  oncogenesis.  Nature.  2017  June  29;546(7660):671-5.  Published. 
Federal  support  was  acknowledged. 

Presentations: 

Bose  R,  Karthaus  W,  Armenia  J,  Abida  W,  Iaquinta  P,  Wongvipat  J,  Wasmuth  E,  Shah  N,  Doran 
M,  Wang  P,  Sullivan  P,  Patruno  A,  International  SU2C/PCF  Consortium,  Zheng  D,  Schultz  N, 
Sawyers  CS.  Loss  of  function  mutations  in  ETS2  repressor  factor,  ERF,  reveal  a  balance  between 
positive  and  negative  ETS  factors  controlling  prostate  oncogenesis.  Prostate  Cancer  Foundation 
Annual  Scientific  Retreat  -  Carlsbad  (2016) 


Website(s)  or  other  Internet  site(s) 
Nothing  to  report 


Technologies  or  techniques 
Nothing  to  report 


Inventions,  patent  applications,  and/or  licenses 
Nothing  to  report 


Other  Products 
Nothing  to  report 
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PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 


Individuals 


Name:  Rohit  Bose 

Project  Role:  PI 

ORCID  ID:  0000-0002-6785-0697 

Nearest  person  month  worked:  1 0.8  months 

Contribution  to  Project:  PI,  Designed,  performed  and  analyzed  experiments. 


Change  in  Other  Active  Support 
Nothing  to  report 


Other  Organizations 
Nothing  to  report 


SPECIAL  REPORTING  REQUIREMENTS 

None 
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LETTER 


doi:10.1038/nature22820 


ERF  mutations  reveal  a  balance  of  ETS  factors 
controlling  prostate  oncogenesis 

Rohit  Bose1,2,  Wouter  R.  Karthaus1,  Joshua  Armenia1,3,  Wassim  Abida2,  Phillip  J.  Iaquinta1,  Zeda  Zhang1,4,  John  Wongvipat1, 
Elizabeth  V.  Wasmuth1,  Neel  Shah1,4,  Patrick  S.  Sullivan1,  Michael  G.  Doran1,  Ping  Wang5,  Anna  Patruno2,  Yilin  Zhao5, 
International  SU2C/PCF  Prostate  Cancer  Dream  Team* *,  Deyou  Zheng5,6,7,  Nikolaus  Schultz1,3, &  Charles  L.  Sawyers1,8 


Half  of  all  prostate  cancers  are  caused  by  the  TMPRSS2-ERG  gene- 
fusion,  which  enables  androgens  to  drive  expression  of  the  normally 
silent  E26  transformation-specific  (ETS)  transcription  factor  ERG  in 
prostate  cells1,2.  Recent  genomic  landscape  studies  of  such  cancers3-8 
have  reported  recurrent  point  mutations  and  focal  deletions  of 
another  ETS  member,  the  ETS2  repressor  factor  ERE9.  Here  we  show 
these  ERE  mutations  cause  decreased  protein  stability  and  mostly 
occur  in  tumours  without  ERG  upregulation.  ERE  loss  recapitulates 
the  morphological  and  phenotypic  features  of  ERG  gain  in  normal 
mouse  prostate  cells,  including  expansion  of  the  androgen  receptor 
transcriptional  repertoire,  and  ERE  has  tumour  suppressor  activity 
in  the  same  genetic  background  of  Pten  loss  that  yields  oncogenic 
activity  by  ERG.  In  the  more  common  scenario  of  ERG  upregulation, 
chromatin  immunoprecipitation  followed  by  sequencing  indicates 
that  ERG  inhibits  the  ability  of  ERF  to  bind  DNA  at  consensus  ETS 
sites  both  in  normal  and  in  cancerous  prostate  cells.  Consistent  with 
a  competition  model,  ERF  overexpression  blocks  ERG-dependent 
tumour  growth,  and  ERF  loss  rescues  TMPRSS2-ERG- positive 
prostate  cancer  cells  from  ERG  dependency.  Collectively,  these 
data  provide  evidence  that  the  oncogenicity  of  ERG  is  mediated, 
in  part,  by  competition  with  ERF  and  they  raise  the  larger  question 
of  whether  other  gain-of-function  oncogenic  transcription  factors 
might  also  inactivate  endogenous  tumour  suppressors. 

Recent  exome  sequencing  revealed  that  3%  of  patients  in  the  SU2C- 
294  metastatic  prostate  cancer  cohort3  (Extended  Data  Fig.  la)  have 
somatic  point  mutations,  but  not  gene  fusions,  involving  the  ETS 
member  ERF.  ERF  was  one  of  a  small  number  of  genes  whose  muta¬ 
tion  frequency  and  predicted  functional  impact  reached  significance  by 
the  MutSig  algorithm10.  Loss-of-function  ERF  germline  mutations  and 
lower  ERF  expression  have  been  previously  implicated  in  the  disease 
complex  craniosynostosis11.  Notably,  the  DNA-binding  ETS  domain  of 
ERF  is  most  similar  to  the  ERG  subfamily12.  However,  unlike  oncogenic 
ETS  factors12,  ERF  possesses  a  transferable  carboxy- terminal  domain 
that  mediates  transcriptional  repression9. 

We  queried  additional  prostate  cancer  genome  cohorts  ( n  =  930 
patients)4-8  (Extended  Data  Fig.  la)  and  found  further  evidence  of 
ERF  mutations  (Fig.  la)  in  1-3%  of  patients.  The  mutations  include 
the  specific  K401fs  and  G299fs  loss-of-function  truncations  also  found 
in  craniosynostosis  families11,  as  well  as  similar  missense  mutations  in 
the  ETS  domain.  Mapping  of  these  ETS  missense  mutations  onto  the 
known  crystal  structure  of  ERG13  revealed  that  the  altered  residues 
are  located  within  conserved  helices  (Extended  Data  Fig.  lb),  with 
four  of  the  five  conserved  mutations  predicted  to  be  destabilizing14. 
Expression  of  complementary  DNAs  (cDNAs)  containing  the  ETS 
missense  or  truncating  mutations  in  LNCaP  prostate  cancer  cells  led 


to  reduced  ERF  steady- state  levels  relative  to  wild- type  ERF  despite 
robust  messenger  RNA  (mRNA)  expression  (Extended  Data  Fig.  lc,  d), 
in  agreement  with  the  craniosynostosis  studies11.  In  addition,  only  full- 
length  ERF  was  detected  following  immunoprecipitation  in  prostate 
cancer  cells  possessing  an  endogenous  heterozygous  K401fs  allele5 
(Extended  Data  Fig.  le).  Consistent  with  destabilization,  we  were 
unable  to  isolate  mutant  ERF  ETS  domains  through  recombinant 
expression  in  bacteria,  whereas  appreciable  amounts  of  wild-type  ERF 
ETS  were  obtained  (Extended  Data  Fig.  If). 

In  addition  to  destabilizing  mutations,  we  observed  remarkably 
narrow  deletions  of  the  ERF  locus  within  the  TCGA-333  primary 
tumour  cohort4  (Fig.  lb).  Median  ERF  expression  in  those  tumours 
containing  focal  hemizygous  deletions  is  lower  than  in  normal 
prostate  and  diploid  ERF  tumours  (Extended  Data  Fig.  2a,  P  =  0.019). 
Intriguingly,  tumours  with  either  ERF  mutations  or  focal  deletions 
are  mostly  exclusive  to  TMPRSS2-ERG-negative  tumours  (Fig.  lc, 
P  =  0.022).  Metastatic  tumours  in  the  SU2C-294  cohort3  containing 
ERF  point  mutations  are  also  mostly  exclusive  to  those  without 
upregulated  ERG,  but  this  distinction  does  not  reach  statistical 
significance  (Extended  Data  Fig.  2b,  P  =  0.066). 

The  near  complete  lack  of  ERF  mutations  in  tumours  with 
TMPRSS2-ERG  fusions  led  us  to  investigate  whether  ERF  loss  reca¬ 
pitulates  the  phenotype  of  ERG  gain.  ERG  is  not  expressed  in  benign 
prostate  epithelium  and  displays  potent  oncogenicity  in  a  mouse 
Pten background2,15.  On  the  other  hand,  ERF  is  endogenously 
expressed  in  normal  prostate  (Extended  Data  Fig.  2a).  To  determine 
whether  its  loss  yields  a  phenotype  similar  to  ERG  gain,  we  infected 
prostate  organoids16  derived  from  Pten+,+  and  Pten~'~  mice15  with 
short  hairpin  RNA  (shRNA)  targeting  mouse  Erf( shErf_m)  (Extended 
Data  Fig.  3a).  The  Pten+I+  shErf_m  organoids  acquired  morphological 
characteristics  of  ERG  overexpression17:  they  formed  single-cell  luminal 
structures  lacking  basal  cells  (Fig.  2a),  and  profiling  by  RNA  sequencing 
(RNA-seq)  demonstrated  profound  loss  of  expression  of  most  basal 
signature  genes18  (Fig.  2b  and  Extended  Data  Fig.  2b-d).  Consistent 
with  these  changes17,18,  RNA-seq  profiling  revealed  significant  enrich¬ 
ment  for  genes  whose  expression  is  upregulated  by  androgen  in  human 
prostate  cancer  cells  (Extended  Data  Fig.  2c,  e).  The  Pten~'~  organoids16 
infected  with  shErf_m  also  demonstrated  a  marked  luminal  shift 
(Fig.  2a)  and  were  able  to  form  tumours  when  grafted  back  into  mice, 
recapitulating  the  phenotype  of  ERG  overexpression15  (Fig.  2c). 

To  assess  the  impact  of  ERF  loss  in  the  half  of  human  prostate  cancers 
that  lack  ERG  expression  similarly  to  normal  prostate,  we  infected  a 
human  cancer-derived  cell  line  lacking  the  TMPRSS2-ERG  gene  fusion, 
CWR22Pc,  with  shRNA  targeting  ERF  (Extended  Data  Fig.  4a)  and 
analysed  its  androgen  transcriptome.  We  observed  both  an  increase  in 
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Figure  1  |  Recurrent  ERF  loss-of-function  mutations  and 
focal  deletions  are  found  in  prostate  cancer  and  are  mostly 
exclusive  to  tumours  without  TMPRSS2-ERG.  a,  ERF 

point  mutations  identified  in  all  prostate  cancer  cohorts3-8. 

*  Mutations  shared  with  patients  having  craniosynostosis. 
Green:  ETS  missense  mutations  of  residues  conserved  with 
ERG.  Red:  frameshift  (fs)/splice-site  (sp)  mutations  upstream 
of  repressor  domain,  b,  ERF  copy  number  deletions  in 
TCGA-333  human  primary  prostate  cancer  cohort4, 
c,  cBio  Oncoprint  of  patients  with  primary  prostate  cancer 
( n  =  333  patients)4.  Each  column  corresponds  to  a  unique 
patient’s  tumour  profile.  P  value  calculated  by  Fishers  exact 
two-tailed  test. 


the  number  of  differentially  expressed  androgen  receptor  target  genes 
and  in  the  magnitude  of  the  expression  changes  (Extended  Data  Fig.  4b), 
despite  no  alteration  of  androgen  receptor  mRNA  or  protein  levels 
(Extended  Data  Fig.  4a,  c).  Next,  we  interrogated  the  mRNA  expres¬ 
sion  profiles  of  the  primary  and  metastatic  human  prostate  cohorts 
TCGA-333  and  SU2C-150  (refs  3, 4),  respectively.  In  agreement  with 
our  functional  studies,  ERF  mRNA  levels  are  inversely  correlated  with 
two  androgen  transcriptional  activity  signatures19,20,  both  in  normal 
human  prostate  and  in  all  primary  tumour  subtypes  (Fig.  2d  and 


Extended  Data  Fig.  5a).  This  reciprocal  association  is  also  observed  in 
metastatic  cancers  if  the  analysis  is  limited  to  tumours  without  ampli¬ 
fication  of  or  mutations  in  the  AR  gene,  which  encodes  the  androgen 
receptor  (Extended  Data  Fig.  5b).  The  fact  that  reduced  ERF  expres¬ 
sion  enhances  androgen  receptor  transcriptional  output  even  in  the 
absence  of  mutation  or  deletion,  raises  the  possibility  that  ERF  may 
have  a  broader  role  in  prostate  oncogenesis. 

The  inverse  correlation  between  ERF  mRNA  level  and  the  andro¬ 
gen  receptor  signature  is  also  observed  in  the  half  of  prostate  cancers 
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Figure  2  |  ERF  is  a  tumour  suppressor,  a  negative 
regulator  of  androgen  signalling,  and  its  loss  phenocopies 
TMPRSS2-ERG  gain,  a,  Immunohistochemistry  of  mouse 
prostate  organoids  infected  with  non-targeting  shRNA 
(shNT)  or  targeting  ERF  (shErf_m).  b,  Basal  signature18 
applied  to  Pten+I+  organoid  RNA-seq  (n  =  2  biological 
replicates),  c,  Tumour  volumes  of  organoid  grafts.  Data  are 
median  ±  interquartile  range;  n  =  10  tumours  per  condition, 
d,  Human  TCGA  primary  prostate  cancer  (PCa)  cohort4 
(n  =  333  patients)  interrogated  for  ERF  mRNA  level  and  the 
Hieronymus  androgen  transcriptional  activity  signature4,19,20. 
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Figure  3  |  ERF  and  TMPRSS2-ERG  have  opposing  effects  on  the 
androgen  receptor  transcriptome  and  cistrome.  a,  Androgen-regulated 
genes  (as  in  Extended  Data  Fig.  4b)  analysed  by  magnitude  of  change 
in  VCaP  cells  infected  with  a  doxycycline-inducible  shRNA  targeting 
ERG  (with  doxycycline,  ERG-low;  without  doxycycline,  ERG-high)  or  a 
constitutive  shRNA  targeting  ERF  (shERF_2);  n  =  3  biological  replicates. 
DHT  downreg.,  DHT  downregulated.  b,  ERF  ChIP-seq  in  VCaP  cells 
(n  =  2  biological  replicates:  Rl,  R2)  analysed  by  (1)  peak  overlap,  (2)  motif 
analysis,  (3)  overlap  with  ERG  and  androgen  receptor  ChIP-seq,  and 
(4)  example  ChIP-seq  signals,  c,  ERF  ChIP-seq  in  normal  prostate 
organoids  infected  with  either  a  Flag-ERG  lentivirus  or  Flag  alone, 
both  treated  with  doxycycline  (n  =  2  biological  replicates:  Rl,  R2)  and 
analysed  as  in  b. 


15  20 

Time  after  injection  (weeks) 


Normal  prostate  Prostate  cancer 

Figure  4  |  TMPRSS2-ERG  activity  is  mediated,  in  part,  by  inactivation 
of  ERF  function,  a,  Tumour  volumes  of  grafts  derived  from 
Pten~l~;R26ERG/ERG  organoids  infected  with  CRISPR-Cas9  targeting 
ERF  (sgERF).  Data  are  median  ±  interquartile  range;  n  =  10  tumours  per 
condition,  P  value  via  Mann- Whitney  exact  two-tailed  test,  b,  Similar  to  a 
but  infected  instead  with  doxycycline-inducible  Flag-ERF.  c,  Cell  viability 
assay  in  pooled  VCaP  cells  infected  first  with  shRNA  targeting  ERF 
(shERF_2)  followed  by  sgERF.  Data  are  mean  =b  s.e.m.;  n  =  3  biological 
replicates,  d,  Androgen-regulated  genes  (as  in  Extended  Data  Fig.  4b)  in 
ERG-low  VCaPs  and  shERF_2  analysed  by  magnitude  of  change;  n  =  3 
biological  replicates.  DHT  downreg.,  DHT  downregulated. 
e,  ERF  inactivation  model  for  genes  with  androgen  receptor  and 
ETS  binding  sites. 


that  possess  the  TMPRSS2-ERG  fusion  (Fig.  2d),  suggesting  that  ERF 
may  have  an  androgen  receptor  repressive  function  in  this  subtype 
as  well.  To  investigate  this  possibility,  the  expression  of  ERG  or  ERF 
was  separately  inhibited  via  shRNA  (Extended  Data  Fig.  6a,  b)  in  the 
ERG-fusion-positive  VCaP  cell  line  and  the  androgen  transcriptome 
analysed  as  before.  Consistent  with  earlier  work15,  inhibition  of  ERG 
expression  (‘ERG-low’)  resulted  in  a  contracted  androgen  transcrip¬ 
tome  compared  with  the  wild-type  state  (‘ERG-high’).  Conversely,  ERF 
inhibition  increased  the  change  in  expression  of  androgen  receptor 
target  genes  and  doubled  the  size  of  the  androgen  transcriptome  (Fig.  3a 
and  Extended  Data  Fig.  7).  ERF  and  ERG  knockdown  had  no  effect 
on  each  other’s  expression,  on  androgen  receptor  protein  levels,  or  on 
androgen  receptor  subcellular  localization  (Extended  Data  Fig.  6c,  d). 

Given  the  similarity  of  their  ETS  domains12  (Extended  Data  Fig.  lb), 
we  postulated  that  the  opposing  effects  of  ERF  and  ERG  on  androgen 
signalling  could  be  explained  by  competition  for  androgen -receptor- 
associated  ETS  binding  sites,  which  we  investigated  by  ERF  chromatin 


immunoprecipitation  followed  by  sequencing  (ChIP-seq)  in  the  ERG- 
high  and  ERG-low  states  (Fig.  3b  and  Extended  Data  Fig.  8).  De  novo 
motif  analysis  identified  the  canonical  ETS  motif  as  the  primary  ERF 
binding  site,  with  2,793  binding  sites  in  the  ERG-high  condition  ((1) 
and  (2)  in  Fig.  3b).  Remarkably,  an  additional  26,714  ERF  binding  sites 
were  observed  in  the  ERG-low  state,  76%  of  which  were  bound  by  ERG 
before  ERG  knockdown  ((1)  and  (3)  in  Fig.  3b,  and  Extended  Data 
Fig.  8d).  Furthermore,  the  ChIP-seq  signal  intensity  of  the  smaller 
number  of  ERF  peaks  observed  in  the  ERG-high  state  was  increased 
in  the  ERG-low  state  in  almost  all  cases  (Extended  Data  Fig.  8c,  e). 
Finally,  ERF  peaks  in  the  ERG-low  state  largely  overlap  with  androgen 
receptor  binding  peaks  ((3)  in  Fig.  3b),  as  previously  reported  for 
ERG15,  and  are  exemplified  by  PLEKHD1  and  SCD20  ((4)  in  Fig.  3b, 
and  Extended  Data  Fig.  8f).  In  the  ERG-high  state,  both  genes  have 
prominent  androgen-receptor-  and  ERG-associated  binding  but  only 
limited  ERF  binding.  However,  in  the  ERG-low  state,  ERF  binding  is 
substantially  increased  in  parallel  with  a  decrease  in  their  androgen- 
induced  expression  (Fig.  3a). 
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Given  the  role  that  ERF  also  plays  in  benign  prostate  cells,  we  per¬ 
formed  ChIP-seq  to  confirm  that  ERF  binds  ETS  sites  in  the  normal 
prostate  organoids  (Fig.  3c).  Moreover,  transient  overexpression  of 
ERG  led  to  a  significant  decrease  of  ERF  chromatin  occupancy  (Fig.  3c 
and  Extended  Data  Fig.  9),  consistent  with  the  competition  for  binding 
also  seen  in  VCaP  cancer  cells.  We  explored  the  overlap  between 
androgen  receptor  and  ERF  sites  in  normal  prostate  organoids  by 
androgen  receptor  ChIP-seq  and  found  28%  of  the  ERF  sites  overlap 
with  androgen  receptor  binding  sites  ((3)  in  Fig.  3c).  This  was  lower 
than  observed  in  the  ERG-positive  VCaP  tumour  cells,  consistent 
with  differences  in  androgen  receptor  chromatin  occupancy  between 
normal  prostate  and  cancer21. 

We  next  asked  whether  ERF  could  modulate  oncogenicity  mediated 
by  ERG  expression.  First,  we  used  CRISPR-Cas9  (clustered  regularly 
interspaced  short  palindromic  repeats  and  CRISPR- associated  protein  9) 
to  partially  delete  ERF  (sgErf)  in  pooled  mouse  Pten~/~;R26ERG/ERG  orga¬ 
noids  (Extended  Data  Fig.  10a),  which  require  ERG  to  form  tumours15. 
The  sgErf  organoids  formed  tumours  more  rapidly  than  those  infected 
with  a  non-targeting  CRISPR  (Fig.  4a).  Likewise,  induction  of  human  ERF 
expression  in  these  cells  blocked  androgen- dependent  gene  expression 
(Extended  Data  Fig.  10b,  c)  and  prevented  tumour  formation  (Fig.  4b). 
To  address  the  role  of  ERF  in  transformation  of  human  prostate  cancer 
cells,  we  predicted  that  the  anti-proliferative  effect  of  ERG  knockdown  in 
ERG-positive  cells22  should  be  rescued  by  concurrent  knockdown  of  ERF. 
Consistent  with  previous  findings22,  partial  ERG  deletion  via  CRISPR- 
Cas9  led  to  a  complete  halt  of  proliferation  (Fig.  4c  and  Extended  Data 
Fig.  lOd).  Strikingly,  stable  shRNA  knockdown  of  ERF  before  partial  ERG 
deletion  rescued  cells  from  ERG- dependent  proliferation  and  survival. 
Moreover,  ERF  knockdown  restored  the  constricted  androgen  transcrip - 
tome  conferred  by  ERG  loss,  including  restoration  of  SCD  and  PLEKHD1 
mRNA  upregulation  (Fig.  4d  and  Extended  Data  Fig.  lOe). 

We  propose  that  loss  of  ERF  activity,  either  by  rare  genomic  loss- 
of-function  mutations  or  more  commonly  by  competition  with  the 
TMPRSS2-ERG  oncogenic  gene  product,  leads  to  activation  of  the 
androgen  receptor  pathway  and  prostate  cancer  (Fig.  4e).  The  con¬ 
sequences  of  ERF  loss  and  ERG  gain  have  many  common  functional 
consequences,  including  regulation  of  luminal  morphology,  expanded 
activity  of  androgen  receptor,  and  an  ability  to  form  tumours  in  a 
PtenT1-  background.  Such  similarities  are  consistent  with  their  mutual 
exclusivity  and  the  ability  of  ERF  loss  to  rescue  ERG-positive  prostate 
cancer  cells  from  ERG  dependency.  On  the  other  hand,  the  fact  that 
ERG  translocations  are  more  common  than  ERE  mutant  tumours 
(46%  versus  4%,  respectively,  in  the  TCGA-333  primary  prostate 
cancer  cohort)  raises  the  possibility  that  ERG  has  additional  gain-of- 
function  activities  that  favour  oncogenic  transformation,  such  as  an 
intrinsic  androgen  receptor  reprogramming  activity15,21.  A  related  but 
unresolved  question  is  whether  the  oncogenic  phenotype  of  ERF  loss 
is  simply  due  to  loss  of  repression,  or  to  gain -of- function  conferred 
by  competition  with  an  unknown  endogenous  positive  ETS  factor. 
Understanding  ETS  competition  has  implications  also  for  ETS-driven 
leukaemias  and  sarcomas,  and  perhaps  more  broadly  for  other  onco¬ 
genic  transcription  factors  for  which  corresponding  loss-of-function 
mutations  in  cognate  repressors  might  be  considered. 

Online  Content  Methods,  along  with  any  additional  Extended  Data  display  items  and 
Source  Data,  are  available  in  the  online  version  of  the  paper;  references  unique  to 
these  sections  appear  only  in  the  online  paper. 
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METHODS 

Cell  lines.  Mouse  prostate  organoids  were  isolated  and  grown  as  described 
previously16  from  the  prostates  of  normal  (Pten+I+),  PtenT /_,  and  Pten R26erg/erg 
mice15.  Isolation  and  growth  of  MSK-PCa3  cells  were  described  previously5. 
LNCaP  cells  were  obtained  from  American  Type  Culture  Collection  (ATCC)  and 
maintained  as  previously  described23.  CWR22Pc  were  a  gift  from  M.  Nevalainen 
and  maintained  as  described  previously23.  For  CWR22Pc  cells  specifically, 
charcoal-stripped  media  was  used  for  transcriptome  analysis.  VCaP  cells  were 
obtained  from  ATCC  (CRL-2876)  and  maintained  as  previously  described23.  Cells 
were  confirmed  to  be  free  of  mycoplasma  using  a  Lonza  detection  kit  (LT07-318). 
Cell  line  authentication  was  confirmed  by  SNP  fingerprinting  and  hallmark 
gene-fusions/ mutations  were  identified  by  deep  sequencing. 

Prostate  cancer  tumour  profiling.  Profiles  of  the  various  cohorts3-7  outlined  in 
Extended  Data  Fig.  la  can  be  explored  in  the  cBioPortal  for  Cancer  Genomics 
(http://www.cbioportal.org).  The  TCGA  data  can  also  be  accessed  through  the 
Broad  Institute  FireBrowse  portal  (http://firebrowse.org/?cohort=PRAD),  and  the 
ICGC-CRUK  Prostate  Adenocarcinoma  data  can  be  accessed  through  the  ICGC 
Portal  (https://icgc.org/icgc/cgp/70/508/71331).  In  Fig.  la  and  Extended  Data 
Fig.  2b,  we  report  four  previously  unpublished  mutations  in  ERF  from  a  larger 
Stand  Up  to  Cancer3  (SU2C)  cohort  that  will  be  described  in  the  future. 
Organoid  histology  and  xenografts.  Organoids  were  infected  with  the  indicated 
lentivirus,  selected  with  puromycin  and  further  selected  with  FACS.  For 
histology  and  immunohistochemistry,  organoids  were  processed  as  described 
previously5,16.  Immunohistochemistry  was  performed  using  a  Ventana  BenchMark 
ULTRA.  The  anti-CK8/18  antibody  was  purchased  from  Abeam  (ab53280). 
In  vivo  xenograft  experiments  were  performed  as  described  previously23,  using 
7-week-old  male  C.B17  SCID  mice  (Taconic):  one  million  cells  were  injected  into 
the  flank  for  a  total  of  ten  tumours  per  shRNA  and  cell  type.  Once  tumours  were 
palpable,  tumour  volume  was  measured  weekly  using  a  Peira  TM900  system  (Peira, 
Belgium).  The  maximal  tumour  volume  permitted  by  our  Institutional  Animal 
Care  and  Use  Committee  under  protocol  06-07-012  was  2  cm3,  beyond  which 
mice  were  euthanized.  All  animal  experiments  were  performed  in  compliance 
with  the  guidelines  of  the  Research  Animal  Resource  Center  of  Memorial  Sloan 
Kettering  Cancer  Center. 

Inhibition  of  gene  expression,  and  overexpression.  ERF  shRNA  knockdown  experi¬ 
ments  were  performed  by  infection  and  puromycin  selection  of  cells  with  lentivirus 
containing  the  miR-E  based  SGEP  or  SGEN  vectors24  gifted  by  J.  Zuber  containing 
the  following  guide  sequences:  shErf_m:  TTGAACTTGTAGGTGAACCGTT, 
shERF_2:  TTGTTTTGAATACATTCTCCAG,  shERF_l:  TTGAAATTGAACT 
TGTAGGTGA,  shNT  was  previously  described  as  Ren.713  targeting  Renilla 
luciferase24.  ERG  shRNA  knockdown  experiments  were  performed  by  infection 
and  neomycin  selection  of  cells  with  lentivirus  containing  the  Tet-pLKO- neo 
vector  gifted  by  D.  Wiederschain  (Addgene  plasmid  21916)  using  the  previ¬ 
ously  published  ERG  shRNA  and  non-targeting  sequences25.  For  Tet-responsive 
reporters,  either  dimethylsulfoxide  (DMSO)  vehicle  or  doxycycline  (Sigma)  was 
added  at  100  ng  ml-1.  CRISPR-Cas9  experiments  were  performed  by  infection 
and  puromycin  selection  of  cells  with  lentivirus  containing  the  lentiCRISPRv2 
vector  gifted  by  F.  Zhang  (Addgene  plasmid  52961)  containing  the  following  guide 
sequences  chosen  via  the  http://www.genome-engineering.org  website:  sgERG 
(GATAACTCTGCGCTCGTTCG),  sgERFm  (CCTGCCAAGCGATGACGCCC), 
and  previously  described  sgNT26.  Overexpression  of  cDNAs  was  achieved  by  the 
constitutive  or  Tet-inducible  pLV-based  lentiviral  expression  system. 
Transcription  analysis.  RNA  was  extracted  from  cell  lines  using  an  RNeasy 
Kit  (Qiagen).  For  quantitative  PCR  with  reverse  transcription  (RT-qPCR)  and 
RNA-seq  experiments,  cells  were  plated  in  triplicate  per  condition/infected 
construct  at  the  beginning  of  the  assay  (duplicate  for  mouse  organoid  RNA-seq), 
and  thereafter  replicates  were  processed  independently.  Error  bars  for  bar 
graphs  indicate  s.e.m.  for  the  biological  replicates.  DHT  (Sigma)  (or  DMSO 
vehicle)  was  added  at  1  nM,  and  treatments  were  performed  for  16  h.  For 
RT-qPCR,  cDNA  was  generated  with  a  High  Capacity  cDNA  Reverse 
Transcription  Kit  (Applied  Biosystems).  Data  were  quantified  relative  to 
(3-actin  or  GAPDH  expression,  and  relative  expression  was  plotted.  Primers  for 
human  and  mouse  ERF  were  purchased  from  Qiagen.  Other  qPCR  primers  were 
as  follows:  androgen  receptor  (F:  CCATCTTGTCGTCAATGTTATGAAGC, 
R:  AGCTTCTGGGTTGTCTCCTCAGTGG),  ERG  (F:  C  A  A  A  ACT  CT  CCA 
CGGTTAATGC,  R:  ACCGGTCCAGGCTGATCT). 

For  RNA-seq,  library  preparation  and  sequencing  were  performed  by  the 
Memorial  Sloan  Kettering  Cancer  Center  Integrated  Genomics  Operation  Core 
using  Illumina  HiSeq  with  50  base  pair  (bp)  paired-end  reads,  with  approximately 
30  million  reads  generated  for  each  sample.  The  output  data  (FASTQ  files)  were 
mapped  to  the  target  genome  (UCSC  HG19  or  UCSC  MM10)  using  the  rnaStar 
aligner,  which  mapped  reads  genomically  and  resolved  reads  across  splice  junctions. 


The  2-pass  mapping  method  was  used  in  which  the  reads  were  mapped  twice.  The 
first  mapping  pass  used  a  list  of  known  annotated  junctions  from  Ensembl.  Novel 
junctions  found  in  the  first  pass  were  then  added  to  the  known  junctions  and  a 
second  mapping  pass  was  done  (on  the  second  pass  the  RemoveNoncanonical 
flag  was  used).  After  mapping,  we  post-processed  the  output  SAM  files  using  the 
PICARD  tools  to  add  read  groups,  with  AddOrReplaceReadGroups,  which  sorted 
the  file  and  converted  it  to  the  compressed  BAM  format.  We  then  computed  the 
expression  count  matrix  from  the  mapped  reads  using  HTSeq  (http://www-huber. 
embl.de/users/anders/HTSeq).  The  raw  count  matrix  generated  by  HTSeq  was 
then  processed  using  the  R/Bio conductor  package  DESeq  (http://www-huber.embl. 
de/users/anders/DESeq/),  which  was  used  both  to  normalize  the  full  dataset  and  to 
analyse  differential  expression  between  sample  groups.  Androgen-regulated  genes 
were  defined  as  a  twofold  difference,  false  discovery  rate  (FDR)  <  0.05  with  1  nM 
DHT  treatment  for  16  h  DHT.  For  GSEA,  statistical  analysis  was  performed  with 
publicly  available  software  from  the  Broad  Institute  (http://www.broadinstitute. 
org/gsea/index.jsp). 

ChIP.  Chromatin  processing,  as  well  as  anti-androgen  receptor,  anti-ERG,  and 
immunoglobulin- G  ChIP  were  described  previously15.  ERF  ChIP  was  per¬ 
formed  with  anti-ERF  antibody  (Pierce  PA5-30237).  For  ChIP-qPCR,  the  ETS2 
promoter  primers  (forward:  TTACTTCCTCCAGAGACTGACGA;  reverse: 
CGCCGGCCAGAGACGAT)  were  used.  The  PSA  upstream  sequence  lacking 
the  ERF  motif  or  androgen  receptor  binding  motif  was  described  previously27. 
For  ChIP-seq,  library  preparation  and  RNA-seq  were  performed  by  the 
Memorial  Sloan  Kettering  Cancer  Center  Integrated  Genomics  Operation  Core 
using  Illumina  HiSeq  with  50-bp  paired-end  reads.  The  reads  were  aligned  to 
the  human  genome  (hgl9,  build  37)  or  the  mouse  genome  (mmlO,  build  38) 
using  the  program  BWA  (default  parameters)  within  the  PEMapper.  The  soft¬ 
ware  MACS2  (ref.  28)  (-q  0.1)  was  used  for  peak  identification  with  data  from 
ChIP  input  DNAs  as  controls.  Peaks  of  sizes  >100  bp  and  with  at  least  one  base 
pair  covered  by  >  15  reads  were  selected  as  the  final  peaks.  Peaks  from  different 
conditions  were  merged  to  obtain  non-overlapping  genomic  regions,  which  were 
then  used  to  determine  conditional  specific  binding.  Overlapped  peaks  were 
defined  as  those  sharing  at  least  one  base  pair.  Two  replicates  were  performed 
per  condition,  and  peaks  called  in  both  replicates  were  used  as  the  final  peaks  for 
each  condition.  To  generate  heat  maps  depicting  ERF  ChIP-seq  read  density  in 
±2  kilobase  regions  of  the  ERF  peak  summits,  the  same  number  of  ChIP-seq  reads 
from  different  conditions  were  loaded  into  the  software  seqMINER29,  and  the 
resulting  read  density  matrices  were  sorted  by  the  read  densities  in  the  ERG-low 
condition  in  VCaP  cells  or  wild-type  condition  for  the  mouse  prostate  organoids, 
before  colouring.  The  criteria  for  assigning  peaks  to  genes  have  been  described 
previously30,31.  The  MEME-ChIP  software32  was  applied  to  300-bp  sequences 
around  the  peak  summits  for  motif  discovery. 

Immunoprecipitation  and  western  blot  analysis.  Protein  was  extracted  from  cell 
lines  using  M-PER  Reagent  (Thermo  Scientific)  and  quantified  by  BCA  Protein 
Assay  (ThermoFisher  Scientific).  Nuclear/cytoplasmic  fractionation  was  achieved 
with  an  NE-PER  kit  (ThermoFisher  Scientific).  ERF  immunoprecipitation  was 
performed  with  anti-ERF  antibody  (Pierce  PA5-30237).  Western  blots  were  imaged 
using  the  fluorometric-based  LiCor  system,  using  primary  antibodies  against 
ERF  (Santa  Cruz  sc-15435  or  Abeam  ab61108)  and  PTEN  (Cell  Signaling  9188). 
Western  blotting  for  androgen  receptor,  ERG,  GAPDH,  and  actin  was  described 
previously15. 

In  vitro  growth  assay.  Cells  were  plated  in  triplicate  and  assayed  at  the  time  points 
indicated  using  CellTiter-Glo  (Promega).  Viability  was  plotted  normalized  to  day 
1.  Error  bars  indicate  s.e.m. 

Statistics.  Mutual  exclusivity  in  the  cBio  Oncoprints  was  calculated  using  a  two- 
tailed  Fisher’s  exact  test.  For  RT-qPCR  and  RNA-seq  experiments,  cells  were 
plated  in  triplicate  per  condition/infected  construct  at  the  beginning  of  the  assay 
(duplicate  for  mouse  organoid  RNA-seq),  and  thereafter  replicates  were  processed 
independently.  RT-qPCR  bar  graphs  are  plotted  as  mean  =b  s.e.m.  For  RNA-seq, 
differentially  expressed  genes  were  defined  as  a  twofold  difference,  FDR  <  0.05  of 
DESeq- normalized  expression.  For  GSEA,  statistical  analysis  was  performed  with 
publicly  available  software  from  the  Broad  Institute  (http://www.broadinstitute. 
org/gsea/index.jsp).  For  the  correlation  of  the  androgen  transcriptional  signatures 
with  ERF  mRNA  in  human  tumours,  statistical  significance  of  correlations  was 
calculated  by  Spearmans  test.  Mouse  graft  experiments  consisted  of  ten  tumours 
per  condition,  and  a  Mann- Whitney  test  was  used  in  Fig.  4  to  analyse  differences 
in  tumour  volume.  The  sample  size  estimate  was  based  on  our  experience  with 
previous  experiments5,15,23.  No  formal  randomization  process  was  used  to  assign 
mice  to  a  given  organoid  injection,  and  experimenters  were  not  blinded. 

Data  availability.  All  data  are  available  from  the  authors  upon  reasonable  request. 
The  described  RNA-seq  and  ChIP-seq  data  have  been  deposited  in  the  Gene 
Expression  Omnibus  under  accession  number  GSE83653. 
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Cohorts 

Source  of  specimen 

Clinical  State 

Cohort 

size 

#  of  patients  with 

ERF  point  mutations 

ERF  Point 
Mutation  Rate 

Reference 

Notes 

TCGA-333 

Primary  site  (prostate) 

Localized, 
treatment  naive 

333 

1 

0.3% 

Cancer  Genome  Atlas  Research  Network. 

The  Molecular  Taxonomy  of  Pnmary  Prostate 
Cancer.  Cell  163, 1011-1025  (2015). 

TCGA-498 

Primary  site 

Localized, 
treatment  naive 

498 

4 

0.8% 

Broac  Institute  FireBrowse  portal: 
(http:/ffirebrowse.org/?cohort=PRAD) 

Includes 

TCGA-333 

CRUK 

Primary  site 

Mixed 

108 

1 

0.9% 

ICGC  Portal:  httos://icac.ora/icac/CQD/ 
70/508/71331 

SU2C-150 

Metastatic  site 

Castration-resistant, 
metastatic  prostate 
cancer 

150 

4 

2.7% 

Robinson.  O.  et  al.  Integrative  clinical 
genomics  of  advanced  prostate  cancer.  Cell 
161,  1215-1228  (2015). 

SU2C-294 

Metastatic  site 

Castration-resistant, 
metastatic  prostate 
cancer 

294 

8 

2.7% 

Pending  publication  from  International  SU2C/ 
PCF  Consortium. 

includes 

SU2C-150 

FHCRC 

Primary  or  metastatic  site 

Castration-resistant, 
metastatic  prostate 
cancer 

54 

2 

3.7% 

Kumar.  A.  et  al.  Substantial  interindividual  and 
limited  intraindividual  genomic  diversity 
among  tumors  from  men  with  metastatic 
prostate  cancer.  Nature  Medicine  1-13 
(2016). 

TrentcVComeliy 
Broad  2016 

Metastatic  site 

Castration-resistant, 
metastatic  prostate 
cancer 

77 

1 

1.3% 

Beltran.  H.  et  al.  Divergent  clonal  evolution  of 
castration-resistant  neuroendocrine  prostate 
cancer.  Nature  Medicine  22.  298-305  (2016). 

MSK  Organoids 

Local  recurrence  or  metastatic 
site  or  circulating  tumor  cell 

Castration-resistant, 
metastatic  prostate 
cancer 

7 

1 

14.3% 

Gao,  D.  et  al.  Organoid  Cultures  Derived  from 
Patients  with  Advanced  Prostate  Cancer.  Cell 
(2014). 
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Extended  Data  Figure  1  |  Recurrent  ERF  loss-of-function  mutations 
are  found  in  prostate  cancer  and  are  destabilizing.  Related  to  Fig.  1. 
a,  Description  of  the  prostate  cancer  patient  cohorts3-8,  b,  ERG  ETS 
domain  (Protein  Data  Bank  accession  number  4IRI)  illustrating  sequence 
conservation  with  ERF.  ERF  missense  mutations  indicated  by  sticks, 
c,  Expression  of  ERF  mutant  proteins  in  LNCaP  cells.  For  gel  source 


data,  see  Supplementary  Fig.  1.  d,  ERF  mutant  protein  densitometry 
from  immunoblot  in  Extended  Data  Fig.  lc  compared  with  mRNA 
RT-qPCR.  e,  ERF  immunoprecipitation  in  VCaP  and  MSK-PCa3 
cells,  f,  Left,  bacterially  expressed  ETS  domain  (hERF20-112)  with  ERF 
tumour  mutations  before  (— )  or  after  (+)  Ulpl  protease  cleavage.  Right, 
densitometry  of  in  vitro  ETS  domains. 


©  2017  Macmillan  Publishers  Limited,  part  of  Springer  Nature.  All  rights  reserved. 


LETTER 


RESEARCH 


a 


ERF  deletions:  none  none  shallow  shallow  deep 
non-focal  focal 

nr  - 

normal  primary 

prostate  prostate  cancer 


b 


p =0.0656 

_ 1  J_ 

CD/*  ««•/  .liiiiiUiLihrmikikiLiimikuiLitm  mimiiitihiimiiihiMiiiiiiiiiiUiiiuikiiiiihiiiiiiiiihhiimtu 


ERF  3%l 
(pt  muts 
only) 


missonse  ■  truncating  )  ERG-high 

mutation  mutation  (>50th  percentile  expression) 


Extended  Data  Figure  2  |  Recurrent  ERF  loss-of-function  mutations 
and  focal  deletions  are  found  in  prostate  cancer  and  are  mostly 
exclusive  to  tumours  without  TMPRSS2-ERG.  Related  to  Fig.  1.  a,  ERF 
expression  in  TCGA-333  cohort4  (n  =  333  patients)  segregated  by  copy 
number  loss.  Data  are  Tukey  box-and-whisker  plots;  P  value  calculated  by 


two-tailed  t-test  of  log2(RSEM  values),  b,  cBio  Oncoprint  of  the  SU2C-294 
(n  =  294)  metastatic  prostate  cancer  cohort3.  Unlike  Fig.  lc,  focal  ERF 
deletions  could  not  be  ascertained  because  of  differences  between  TCGA 
and  SU2C  copy  number  data3,4.  P  value  calculated  by  Fishers  exact  two- 
tailed  test. 
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Extended  Data  Figure  3  |  ERF  is  a  negative  regulator  of  androgen 
signalling.  Related  to  Fig.  2.  a,  Mouse  prostate  organoids  infected  with 
non-targeting  shRNA  (shNT)  or  targeting  ERF  (shErf_m),  grown  in  three- 
dimensional  culture.  For  RT-qPCR,  n  —  2  biological  replicates.  For  gel 
source  data,  see  Supplementary  Fig.  1.  b,  RT-qPCR  analysis  of  the  Pten+I+ 
organoids;  n  =  2  biological  replicates,  c,  RNA-seq  analysis  of  organoids 


derived  from  Pten+/+  and  Pten~!~  mouse  prostates  infected  with  non¬ 
targeting  shNT  or  shErf_m;  n  =  2  biological  replicates,  d,  Pten+I+  organoid 
RNA-seq  (n  =  2  biological  replicates)  interrogated  by  GSEA  for  expression 
signature  of  Witte  basal  prostate  cancer18,  e,  Same  data  interrogated  by 
GSEA  for  Nelson  androgen  up  expression  signature20. 
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Extended  Data  Figure  4  |  ERF  is  a  negative  regulator  of  androgen 
signalling.  Related  to  Fig.  2.  a,  CWR22Pc  cells  infected  with  shRNA 
targeting  human  ERF  (shERF_l)  or  a  non-targeting  shRNA  (shNT).  For 
gel  source  data,  see  Supplementary  Fig.  1.  b,  Androgen-regulated  genes 
(at  least  a  twofold  change,  FDR  <  0.05  by  RNA-seq  with  1  nM  DHT 


for  16h)  in  CWR22PC  cells  infected  with  either  shERF_l  or  shNT, 
analysed  by  the  number  (left,  Venn  diagram),  the  magnitude  of  expression 
change  (centre,  graph),  and  heat  map  (right);  n  =  3  biological  replicates, 
c,  RT-qPCR  from  CWR22Pc  cells  infected  with  shERF_l  or  shNT. 

Data  are  mean  ±  s.e.m.;  n  =  3  biological  replicates. 
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Extended  Data  Figure  5  |  ERF  is  a  negative  regulator  of  androgen 
signalling.  Related  to  Fig.  2.  a,  Full  version  of  Fig.  2d.  Expression  profiles 
of  the  TCGA-333  primary  prostate  cancer  cohort4  ( n  =  333  patients) 
were  interrogated  for  correlation  between  the  ERF  mRNA  level  and  two 


androgen  transcriptional  activity  signatures4,19.  P  values  were  calculated 
by  the  Spearman  correlation  test,  b,  The  same  analysis  as  a  was  applied  to 
the  SU2C-150  (n  =  150  patients)  metastatic  castration-resistant  prostate 
cancer  cohort3  (mCRPC). 
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Extended  Data  Figure  6  |  ERF  and  ERG  knockdown  do  not  affect 
androgen  receptor  levels  or  its  subcellular  localization.  Related  to 
Fig.  3.  a,  VCaP  cells  infected  with  doxycycline  (dox) -inducible  shRNA 
targeting  ERG  (+dox,  ERG-low;  —dox,  ERG-high).  For  gel  source  data, 
see  Supplementary  Fig.  1.  b,  VCaP  cells  infected  with  shRNA  targeting 


ERF  (shERF_2)  or  a  non-targeting  shRNA  (shNT).  c,  Nuclear/cytoplasmic 
fractionation  of  VCaP  cells  infected  with  shNT  or  shERF_2,  and  the 
doxycycline-inducible  shRNA  targeting  ERG  (shERG).  d,  RT-qPCR  with 
ERG-low  VCaP  cells  compared  with  those  infected  with  shERF_2.  Data  are 
mean  =t  s.e.m.;  n  =  3  biological  replicates. 
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Extended  Data  Figure  7  |  ERF  and  TMPRSS2-ERG  have  opposing 
effects  on  the  androgen  transcriptome.  Related  to  Fig.  3.  a,  Androgen- 
regulated  genes  (at  least  a  twofold  change,  FDR  <  0.05  by  RNA-seq  with 
1  nM  DHT  for  16  h)  in  VCaP  cells  infected  with  a  doxycycline-inducible 
shRNA  targeting  ERG  (with  doxycycline,  ERG-low;  without  doxycycline, 
ERG-high)  or  a  constitutive  shRNA  targeting  ERF  (shERF_2)  and 
analysed  by  number  (top)  and  heat  map  (bottom  left);  n  =  3  biological 


replicates.  Bottom  centre,  RNA-seq  analysis  evaluating  the  effect  of 
dihydrotestosterone  on  ERF  expression.  Bottom  right,  the  effect  of 
doxycycline  alone  on  RNA-seq  differential  expression  analysis, 
b,  RT-qPCR  of  shERF_2-infected  VCaP  cells  treated  =t  DHT.  Data  are 
mean  ±  s.e.m.;  n  =  3  biological  replicates,  c,  Interrogation  of  RNA-seq 
in  shERF_2  VCaP  cells  (n  =  3  biological  replicates)  by  GSEA  for  Nelson 
androgen  up  expression  signature20. 
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Extended  Data  Figure  8  |  ERF  and  TMPRSS2-ERG  have  opposing 
effects  on  the  androgen  receptor  cistrome.  Related  to  Fig.  3. 
a,  ChIP-qPCR  (n  =  2  biological  replicates)  in  VCaP  cells,  amplifying 
either  the  ETS2  promoter  region  that  contains  a  known  ERF  binding 
site9,  or  an  upstream  element  of  PSA27  lacking  the  ERF  binding  motif 
noted  in  (2)  in  Fig.  3b.  b,  The  effect  of  ERF  shRNA  knockdown  on  its 
binding  to  the  ETS2  promoter  as  assessed  by  ChIP-qPCR  (n  =  2  biological 
replicates),  compared  with  its  effect  on  ERF  mRNA  by  RT-qPCR.  Data  are 


mean  ±  s.e.m.;  n  =  3  biological  replicates,  c,  ERF  ChIP-seq  in  ERG-high 
or  ERG-low  (n  =  2  biological  replicates:  Rl,  R2)  analysed  by  heat  maps, 
d,  Comparison  of  ERF  ChIP-seq  peak  numbers,  n  =  2  biological  replicates: 
Rl,  R2.  e,  A  10-Mb  region  illustrating  ChIP-seq  of  ERF  binding  in  ERG- 
high  condition  compared  with  ERG-low.  f,  ChIP-seq  signals  for  the  SCD 
and  PLEKHD1  loci.  In  both  e  and  f,  ChIP-seq  signals  at  they  axis  were 
normalized  by  read  depths. 
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Extended  Data  Figure  9  |  ERG  expression  decreases  the  ERF  cistrome 
in  normal  prostate  organoids.  Related  to  Fig.  3.  a,  Mouse  normal  prostate 
organoids  were  infected  with  a  TetOn  doxycycline-inducible  Flag-ERG  or 
empty  vector  (Flag-alone)  and  treated  with  or  without  doxycycline.  For  gel 
source  data,  see  Supplementary  Fig.  1.  b,  ERF  ChIP-seq  in  normal  prostate 


organoids  infected  with  either  Flag-alone  or  Flag-ERG  lentivirus,  both 
treated  with  doxycycline  {n  =  2  biological  replicates:  Rl,  R2),  and  analysed 
with  heat  maps,  c,  Comparison  of  ERF  ChIP-seq  peak  numbers,  n  =  2 
biological  replicates  (Rl,  R2).  d,  ERF  ChIP-seq  signals  for  the  Rph3al  and 
Sh3bp5  loci,  normalized  by  read  depths. 
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Extended  Data  Figure  10  |  TMPRSS2-ERG  activity  is  mediated,  in 
part,  by  inactivation  of  ERF  function.  Related  to  Fig.  4.  a,  Pooled 
mouse  Pten~'~\  R26erg/erg  organoids  infected  with  CRISPR-Cas9 
targeting  AAVS1  (sgNT)  or  ERF  (sgErf).  For  gel  source  data,  see 
Supplementary  Fig.  1.  b,  R26erg/erg  organoids  infected  instead 

with  a  doxycycline-inducible  Flag-ERF  or  empty  vector  (Flag  alone), 
c,  RT-qPCR  of  mRNA  isolated  from  the  Flag-ERF-infected  organoids  and 
treated  with  or  without  doxycycline,  with  or  without  dihydrotestosterone 
1  nM  for  16  h.  Data  are  mean  zb  s.e.m.;  n  =  3  biological  replicates,  d,  Pooled 


VCaP  cells  were  first  infected  with  non-targeting  shRNA  (shNT)  or  ERF 
(shERF_2),  followed  by  sgNT  or  sgERF  CRISPR-Cas9.  e,  Androgen- 
regulated  genes  (at  least  a  twofold  change,  FDR  <  0.05  by  RNA-seq  with 
1  nM  DHT  for  16h)  in  VCaP  cells  infected  with  a  doxycycline-inducible 
shRNA  targeting  ERG  (with  doxycycline,  ERG-low)  and  either  a 
constitutive  shRNA  targeting  ERF  (shERF_2)  or  a  non-targeting  shRNA 
(shNT).  Data  analysed  by  number  (left,  Venn  diagram)  and  heat  map 
(right);  n  =  3  biological  replicates. 
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